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Abstract 

If sufficiently light dibaryon resonances exist, a Bose condensate of dibaryons 
can occur in nuclear matter before the quark-hadron phase transition. Within a 
relativistic mean-field model we show that heterophase nuclear-dibaryon matter 
is for a wide set of parameters energetically more favorable than normal nuclear 
matter. Production of dibaryons is, however, relatively suppressed as compared to 
estimates based on the model of non-interacting nucleons and dibaryons. 



The possibility of existence of dibaryon resonances was investigated in the last two 
decades experimentally and theoretically. The most promising candidates for experimen- 
tal searches are those dibaryons which have a small width. In 1977 R.Jaffe predicted JI| 
the existence of a loosely bound dihyperon, H, with a mass just below the threshold for 
two-lambda decay. Calculations of the if-particle mass within QC-D-inspired models [2-4] 
showed that the existence of a dihyperon near the AA-threshold is plausible. Dibaryons 
with exotic quantum numbers, which have a small width due to zero coupling to the 
iViV-channel, are of special interest [5-7]. A method for searching narrow, exotic dibaryon 
resonances in the double proton-proton bremsstrahlung reaction is discussed in Ref. [||. 
Data from pion double charge exchange (DCE) reaction on nuclei || [K| exhibit a pe- 
culiar energy dependence at the total pion energy of 190 MeV that can be interpreted 



TTj as evidence for the existence of a narrow d! dibaryon with quantum numbers T = 0, 



J p = 0~. Recent experiments at TRIUMPF (Vancouver) and CELSIUS (Uppsala) seem 
to support the existence of the d! dibaryon [T2|]. 

The properties of nuclear matter with admixture of multiquark clusters are discussed 
in Ref. [T3| . A dibaryon Bose condensate in interiors of neutron stars decreases the 



maximum masses of neutron stars JT4J . In a recent paper JL5] an exactly solvable model 
for a one-dimensional system of fermions interacting through a potential, which leads to 
a resonance in the two-fermion channel, is constructed. The behavior of this system can 
be interpreted in terms of a Bose condensation of the two-fermion resonances. 

There is no dibaryon condensate in ordinary nuclei. From this one can conclude that 
the masses of dibaryons coupled to the iViV-channel should be greater than 

m D > 2fi N = 2(m N + e F ) = 1.96 GeV (1) 

where /i^v is the chemical potential of nucleons and Ef = 40 MeV is the Fermi energy of 
nucleons in nuclei. Here, the ideal gas approximation for nucleons and dibaryons and the 
assumption that the shell model potential for dibaryons is twice as deep as the one for 
nucleons have been used. 

The d! dibaryon is coupled to the NNtt channel only. In the nuclear medium, the 
reaction nd' <-> nnp is possible. In nuclei the equilibrium condition for the chemical 
potentials has the form \i n + mu = 2/i n + fi p . Since fi p ~ /z n , we arrive at the same 
inequality ([!]). 

A Bose condensate of dibaryons can presumably be formed at high densities when 
relativistic effects for nucleons become important. In order to describe such a system, 
we should go beyond non-relativistic many-body theory. The relativistic field-theoretical 
Walecka model |TB[ is known to be very successful in describing properties of infinite 
nuclear matter and of ordinary nuclei throughout the periodic table. In this paper we 
study the influence of narrow dibaryon resonances on nuclear matter in the framework of 
the Walecka model in the mean-field approximation. 

The Lagrangian of the model contains nucleons interacting through u- and a-meson 
exchanges. We add to the Lagrangian dibaryons interacting with nucleons and each other 
through u- and cr-meson exchanges also. Inclusion of dibaryons entails uncertainties 
connected to the lack of reliable information on dibaryon masses and coupling constants. 
However, many conclusions can be drawn on quite general grounds without knowing 
precise values for the newly added parameters. The Lagrangian density is given by 



2^-2 



2 '" a ~ (2) 
~\ F lv + k m l u l + ( 9 f ~ ih^u^ip*^ + ihnujy - (m D + h a a) 2 ip*ip + C c . 

Here, $ is the nucleon field, ui^ and a are fields of the oj- and a-mesons, = dyU^ — d^uiy, 
ip is the dibaryon isoscalar-scalar (or isoscalar-pseudoscalar) field. The values and 



m a are the uj- and cr-meson masses and the values g u , g a , h^, h a are coupling constants 
of the uj- and cr-mesons with nucleons (g) and dibaryons (h). 

The term C c describes conversion of dibaryons into nucleons. The if-particle is coupled 
to the iViV-channel through a double weak decay and C c = 0(G 2 F ). For the nonstrange d\ 
dibaryon |I7j and the d! dibaryon, we neglect possible virtual transitions e.g. to the NNa 
channel. The on-shell couplings for these dibaryons are small. The exotic d% dibaryon 
decays to the iViV7-channel only, and so C c = 0(a). The d' dibaryon decays to the 
NNw channel. Due to Adler's consistency condition Jl8| C c oc d^n . In the mean-field 
approximation <9 M 7r = 0, and the term C c does not modify the mean-field equations. In 
what follows we set C c = 0. The effect of a small term C c reduces to providing a chemical 
equilibrium with respect to transitions between dibaryons and nucleons [Hj . 

The field operators can be expanded in c-numbers and operator parts: uj^ = g^uJc+uJ^, 
a = <j c + a, ip = ip c + 0, and if* = (p* + 0*. The cr-meson mean field determines the effective 
nucleon and dibaryon masses in the medium: m* N = + g a a c and m* D = vtld + h a a c . 
The baryon number current has the form = jjf + 2j® where = 1 J r 7 At \E r and = 
¥>*i dfi f — 2h u uj^ip*ip. The uj- mesons are coupled to the current = g^j^ + h u j® . 

The nucleon vector and scalar densities are defined by the expectation values p^y =< 
^(0)70^(0) > and pns =< ^ (0)\I/(0) >. The scalar density of the dibaryon condensate 
is given by p c DS = \< (f(0) >| 2 . The time evolution of the dibaryon condensate 93-field is 
determined by the dibaryon chemical potential po 



<p e (t) = e^y/fag. (3) 

The vector density of the dibaryon condensate is given by p c DV = 2p* D p c DS where pn = 
P* D + h w uj c . 

The existence of a dibaryon condensate depends on the values of the coupling constants 
of dibaryons with the uj- and a -mesons. The uj- and a- meson coupling constants and 
h a enter the dibaryon-dibaryon Yukawa potential 

l2 -m^r y/1 -m a r 

V(r) = Y~ T-- • W 

Air r 4ir r 

The interaction energy for dibaryons in the condensate is for a constant density distribu- 
tion Pd(x) = pr> equal to 

W= l -J dx 1 rfx 2 p D (x 1 )p D (x 2 )V r (|x 1 - x 2 |) = 2nN D p D (^ I - (5) 

uj a 

where Njj is the total number of dibaryons. A negative W would imply instability of the 
system against compression. The value W is positive and the system is stable for 

h 2 h 2 

>At- (6) 



Airm 2 Anm? 



In a nonrelativistic theory for systems of interacting bosons [^UJ and in the model consid- 
ered, the requirement of stability is equivalent to the requirement of a positive value for 
the square of the sound velocity (a 2 > 0). 

The ii-particle interactions are studied on the basis of the non-relativistic quark cluster 
model [3, 22-24] which is successful in describing the iViV-phase shifts. The calculation of 



the interaction integral (|5|) with the adiabatic if if -potential [23] gives a negative energy, 
so the if-dibaryon condensate is probably unstable against compression. The coupling 
constants of the mesons with the if -particle can be fixed by fitting the depth and the 
position for minimum of the HH -potential to give /i 2 = 603.7 and h 2 = 279.2. 

The coupling constants of the mesons with the d\ and d dibaryons are unknown. The 
uj- and <7-mesons interact with nucleons and pions inside the dibaryon. For dibaryons 
decaying into the iVW-channel, the a — D and uj — D couplings are in nonrelativistic 
approximation two times greater than for nucleons: = 2g UJ and h a = 2g a . The scalar 
charge is, however, suppressed by the Lorentz factor. For standard parameters of the 
Walecka model [fLj, m a = 520 MeV, g* = 190.4, and g 2 a = 109.6, the inequality (6) 
becomes 98.85(^) 2 Ge1/- 2 > 129.0(^) 2 Ge1/~ 2 . With these additive estimates, the in- 
equality (6) is not fulfilled. The exchange current contributions to the meson couplings 
with dibaryons, which violate additivity, are analysed in Ref. [pi]. At present no definite 
conclusions concerning the stability of the d\ and d' dibaryon matter can be drawn. 

The mean-field solutions of the equations of motion corresponding to the Lagrangian 
density (2) are obtained by neglecting the operator parts of the meson fields. For the u- 
and cr-meson mean fields, we get the following expressions 

g^PNY + h^2p* D p c DS 



m 2 , 



(7) 



9apNS + h a 2m* D p c DS 

Substituting expression (|) into the equation of motion for the dibaryon field (p, we get 
H* D = m* D . The nucleon and dibaryon chemical potentials have the form p N = E* F + g^uic 
and pd = fn*o + h^c, where E* F = m*y + k. F is the Fermi energy of nucleons with the 
effective mass m* N . 

The self-consistency condition for the effective nucleon mass can be transformed to a 
form equivalent to that in the standard Walecka model: 

P- 

'? ■ 



m* N = m N - ^p NS , (9) 
mi 



where m N = m N (p c ^ - p c DV )lp c ^ and p%y ax = m N ml/(g a h a ) = 0.1507( 2 ^)/m 3 . 
If the densities pxv and p c DV are fixed, equation (|]) allows to find the effective nucleon 



mass m* N . Solutions to Eq.(^) exist for arbitrary total density ptv an d dibaryon density 
Pdv < P£)™ ax when the value m N is positive. 

It is clear that the fraction of dibaryons should increase when the difference — Pd 
is positive and p c DV = 0. If the difference 2/ijv — Pd is negative and the system consists 
of dibaryons only, production of nucleons is energetically favorable. If the difference 
2fiN — /id = and increases with the dibaryon fraction, small fluctuations take the 
system away from equilibrium. These types of states are unstable. 

There are three stable cases: (i) The homophase nuclear matter: 2p^ — Pd < and 
Pdv = 0; (ii) the homophase dibaryon matter: 2p N — p D > and 2p c DV = p TV ; (iii) 
the heterophase nuclear-dibaryon matter: 2p N — p D = and ^i^zEm |^ < o. Small 
fluctuations around these states lead the system back to the equilibrium points. 

In Fig.l we show the critical density for occurrence of a Bose condensate of non-strange 
dibaryons for = 2g u] as a function of the cx-meson coupling constant with dibaryons. 
The critical density is determined from the equation 2p^ — po — 0. 

In Fig. 2 (a) we show the nucleon effective mass m* N versus the dibaryon fraction 
2p c DV /pTv for = 2g UJ and h a = 1.6^0-. The behavior of 1YI j\j clS cL function of the 
dibaryon fraction does not depend on the mo, since the dibaryon mass does not enter the 
self-consistency condition (|9]) directly. In Fig. 2 (b) we show the difference for the chemical 
potentials versus the dibaryon fraction. 

Mean-field solutions exist at all densities ptv for sufficiently small densities of dibaryons, 
Pdv< Pz)™ ax - This means that we can always investigate the stability of normal nuclear 
matter with respect to dibaryon condensate formation. Small total baryon number densi- 
ties correspond to a stable equilibrium of type (i), at higher densities a stable equilibrium 
(iii) occurs. When the density ptv is high, dibaryon production is energetically favorable. 
The mean-field solutions disappear, however, before the system reaches an equilibrium. 

In Fig. 3 we show the energy per nucleon and the pressure versus the total baryon 
number density for some possible dibaryons. The effect of zero compressibility for het- 
erophase nuclear-dibaryon matter present in the ideal gas approximation reveals itself 
through the softening of the equation of state (EOS). Notice that the pressure of the 
heterophase system obeys the basic inequality ]25| dp/dpxv > of statistical mechanics. 
One can verify that in the model considered, the hydrostatic pressure coincides with the 
thermodynamic pressure. The if -particles are formed at a lower density, since the u — H 
coupling constant h w / {2g UJ ) = 0.89 is relatively small. The energy of the if -particles in the 
positive cu-meson mean field is lower, so the production of the Jf-particles is energetically 
more favourable. 

The qualitative estimates based on a model for non-interacting nucleons and dibaryons 



show that in normal nuclear matter a dibaryon Bose condensate does not exist provided 
the inequality ([I]) is satisfied. A more accurate estimate can be made on the basis of the 
relativistic mean-field model @. From the requirement of absence of a dibaryon Bose 
condensate for ptv < Po = 0.15 fm~ 3 , we get for h u = 2g w a constraint 

m D > 1.89 GeV. (10) 

This constraint is valid provided that the dibaryon matter is stable against compression. 
In such a case, the d\ resonance |ljj with a mass nif) = 1.92 GeV does not affect properties 
of ordinary nuclei. 



Phase transitions of nuclear matter to strange quark matter [26, 27] have been widely 



discussed in the literature (for a recent review see P8]|). Dense nuclear matter with 
a dibaryon Bose condensate can exist as an intermediate state below the quark-gluon 
phase transition. This is the case when dibaryon matter is stable against compression. 
If dibaryon matter is unstable, the creation of dibaryons can be a possible mechanism 
for the phase transition to quark matter. The energetically favourable compression of 
if -matter can lead to the formation of strange matter. It would be interesting to check 
astrophysical data for the presence of a dibaryon condensate in the interiors of massive 
neutron stars as well as possible signatures of their instability caused by dibaryons. 
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FIGURE CAPTIONS 



Fig.l. The critical density for occurrence of a Bose condensate of dibaryons versus 
the (T-meson coupling constant h a for mo = 1.88 GeV (= 2m^; the long dashed curve 
No. 1), 1.96 GeV (the solid curve No. 2), etc. with a step 80 MeV. The results for 
dibaryons di(1920) and d (2060) are also shown (the dashed curves). The dibaryon matter 
is stable against compression when the square of the sound velocity a s is positive. This 
is the case for h a /(2g a ) < 0.8754. The value po — 0.15 fm~ 3 is the saturation density for 
nuclear matter. For 2m^ < mjj < 1.89 GeV, we start at zero density from a heterophase 
nuclear-dibaryon matter. With increasing the density, the matter can be transformed to 
a homophase nuclear matter and then again to a heterophase nuclear-dibaryon matter. 
For vtld > 1.89 GeV, we start at zero density from a homophase nuclear matter which 
converts with increasing the density (at prv > Po for h a f (2g a ) < 0.8754) to a heterophase 
nuclear-dibaryon matter. The occurrence of the if -particles is denoted by the cross. 

Fig. 2. (a) The effective nucleon mass m* N in GeV versus the dibaryon fraction 
2p c DV / ptv m the heterophase nuclear-dibaryon matter. The results do not depend on the 
dibaryon mass, (b) The difference 2p^ — Pd between the two nucleon chemical potentials 
and the dibaryon chemical potential versus the dibaryon fraction 2p c DV / Ptv- The results 
are given for the total baryon densities 1, 2, 3, 4, 5, and 6 times greater than the satura- 
tion density. The normal nuclear matter is stable when 2p^ — prj < and p c DV = 0. An 
intersection of a curve with a negative slope with the horizontal line 2p^ — p^ = in- 
dicates occurrence of a stable equilibrium in the heterophase nuclear-dibaryon matter. 
The results are given for = 1.96 GeV. The dibaryon mass does not enter the self- 
consistency condition (9) and enters linearly in the difference 2p^ — pr>, so the curves for 
other dibaryon masses can be obtained by vertical parallel displacements. The results for 
rriD = 2.06 GeV (cf-dibaryon) can be obtained e.g. by a 100 MeV negative shift, etc. 

Fig. 3. The pressure (left scale) and the energy per baryon (right scale) versus the 
total baryon number density ptv — Pnv + 2p c DV for normal nuclear matter (solid lines) 
and for the heterophase nuclear-dibaryon matter (dashed lines) for the <ii(1920) and 
d'(2060) dibaryons at h u = 2g U) and h a /{2g (T ) = 0.8. The dibaryon condensates occur at 
Ptv I 'po = 2.05 and 3.15 for the d\ and d' . The dibaryon Bose condensation softens the 
EOS for nuclear matter. The if -dibaryons are formed at prv/po — 2.74 when the meson 
coupling constants are fixed by fitting the adiabatic if if-potential |p3|| . In this case the 
if -matter is unstable against compression, providing a transition to strange matter. 



